Chaos synchronization in long-range coupled map lattices 
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We investigate the synchronization phenomenon in coupled chaotic map lattices where the cou- 
plings decay with distance following a power-law. Depending on the lattice size, the coupling 
strength and the range of the interactions, complete chaos synchronization may be attained. The 
synchronization domain in parameter space can be analytically delimited by means of the condition 
of negativity of the largest transversal Lyapunov exponent. Here we analyze in detail the role of 
all the system parameters in the ability of the lattice to achieve complete synchronization, testing 
analytical results with the outcomes of numerical experiments. 
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Coupled map lattices (CMLs), dynamical systems with 
discrete space and time, are being intensively investi- 
gated nowadays as models of spatiotemporal phenomena 
occurring in a wide variety of systems of physical, biolog- 
ical and technical interest [1]. In this letter we will deal 
with the phenomenon of synchronization and, in particu- 
lar, amongst the various kinds of synchronized behavior, 
with the complete synchronization (CS) [2] occurring in 
CMLs with regular long-range interactions. Most work 
done so far on synchronization in CMLs has focused on 
two extreme coupling types: local (nearest-neighbor) [3] 
and global ("mean field") ones [4]. However, non-local 
couplings are relevant to a variety of situations ranging 
from neural networks [5] to physico-chemical reaction sys- 
tems [6]. 

We consider chains of N coupled one-dimensional (ID) 
chaotic maps x i— > f{x) whose evolution is given by [7] 



N' 

.„+i-(l-£)/(a:«) + ^E 



(1) 



where Xn^ represents the state variable for the site i 
{i — 1, 2, N) at time n, e > is the coupling strength, 
a > controls the effective range of the interactions 
and ?7 = 2^^j^r^" is a normalization factor, with 
TV' = (TV - l)/2 for odd TV. The main interest in this 
coupling scheme resides in the fact that it allows to inves- 
tigate the role of the range of the interactions, scanning 
from the local {a oo) to the global (a = 0) cases [8]. 

CS takes place when the dynamical variables that de- 
fine the state of each map adopt the same value for all 
the coupled maps at all times n, i.e, Xn ^ = Xn ^ = . . . = 
Xn = Xn . Depending on the lattice size and on the 
range of the interactions, there may exist an interval of 
values of the coupling strength e, for which such state is 
spontaneously attained, as we have analytically shown in 
previous work [9] . It is our purpose here to scrutinize the 



role of all the system parameters in the ability of the lat- 
tice to synchronize. Analytical results will be compared 
with the outputs of numerical simulations performed for 
diverse chaotic maps. 

Complete synchronization can be characterized by a 
complex order parameter defined, for time n, as i?„ = 

\jjJ2jLi^'^"^"^\ [10]- Typically a time-averaged ampli- 
tude R is computed over a time interval long enough to 
allow the lattice attain the asymptotic state. In the CS 
state, one has ^ = 1 within a small allowed deviation. 

Another diagnostic of complete synchronization can be 
obtained from the Lyapunov spectrum (LS) of the CS 
states. If the chaotic maps are completely synchronized, 
the maximal Lyapunov exponent, in the direction parallel 
to the synchronization manifold (SM), is strictly positive. 
The negativity of the second largest Lyapunov exponent, 
which belongs to the direction transversal to the SM, in- 
dicates the stability of the synchronized state under small 
transversal displacements [11]. 

In our case the lattice dynamics given by Eq. (1) can 
be written as x. 



of the form 



n-l-l 



Si Pijf{xn^)i where F is a matrix 



(l-e)l-F-B 



(2) 



with 1 the N X N identity matrix and B defined by 
Bjk = (1 - Sjk)/rjf,, being rjk = min/g2;|j -k-\-lN\. 

The Lyapunov spectrum is obtained from the dynamics 
of tangent vectors ^, which in turn is obtained by differ- 
entiation of the original evolution equations. In matrix 
form the tangent dynamics reads ^„ = %iCo^ where 
is product of n Jacobian matrices calculated at succes- 
sive points of a given trajectory. If A^^), . . . , are the 

eigenvalues of A = lim (7^ 7^1)2^, the Lyapunov expo- 
rt — ^CxD 

nents are obtained as A*^'^^ = InA*^'^', for k = 1,...,N 
[12]. For the CS states, one arrives at the following ex- 
pression for the Lyapunov spectrum [9] 



1 



Xi'') = Ac/ + In 



1-e + e- 



V 



(3) 



where Aj/ > is the Lyapunov exponent of the uncou- 
pled chaotic map, and b^''^ are the eigenvalues of B that 
can be obtained by Fourier diagonalization and, for odd 

N, read 



N' 



cos{2-Kkm/N) 



l<k< N. 



m" 



(4) 



For even A^, summations run up to A''' = N/2 and half 
of the A'''th term has to be substracted. The maximal 
eigenvalue is fe*^^-* = 77 and the minimal one is Ex- 
cept for the cases k = N' , with even N , and k = N, 
the remaining eigenvalues are two-fold degenerate, being 

^(fc) ^ b(N-k)^ 

In the calculation of Lyapunov exponents A^'^^ notice 
that the parameters that define the particular uncoupled 
map affect only A[/, while the second term in Eq. (3) is 
determined by the particular cyclic dependence on dis- 
tance in the regular coupling scheme (a power law in our 
case). It can be easily verified that, for arbitrary a, the 
CS state lies along the direction of the eigenvector asso- 
ciated to the largest exponent A'^''^^. Therefore, the CS 
state will be transversally stable if the (A^ — 1) remaining 
exponents are negative, that is |1 — e -I- eb'^'^^ /rj] < e~^" , 
\/k ^ N . This is equivalent to requiring that the sec- 
ond largest (or largest transversal) asymptotic exponent, 
denoted by A""-, be negative. This exponent is obtained 
from Eq. (3) with cither k = 1 (hence also k = N —1 due 
to degeneracy) or with k = N' (hence also k = N' + 1 
if N is odd), depending on whether |1 — £ -|- eb^^yrj] is, 
respectively, greater or smaller than |1 — e -I- eb'-^ V*?!- 
The condition A-'* < leads to £c < £ < [9] (see also 
[13]), where 



£e(a,iV,A[/) = (l-e-^") 1 



6(1 



and (5) 



£^(a,iV,Ay) = (l + e-^-)(l-^ . (6) 



V / 

In Fig. 1 we show a variety of critical curves in parame- 
ter space (a, £) obtained for different lattice sizes. Stable 
CS states dwell in the region bounded by the a ~ axis, 
and two curve segments. The critical curves were ob- 
tained analytically from Eqs. (5) (lower curve) and (6) 
(upper curve). The symbols shown stand for the numer- 
ical results determined from the condition ^ = 1 with a 
tolerance of 10~^. Two different values of Xu were consid- 
ered. Numerical results shown in Fig. 1 were computed 
for the piecewise linear (a) Bernoulli shift f{x) = 2x 
(mod 1) (therefore Xu = In 2) and (b) triangular map 
[14] 



fw{x) = 



x/w for < X < w 

{l-x)/{l-w) for u; < x < 1 ^' 



(a) J^Li = ln2 
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FIG. 1. Synchronization diagram in parameter space (a, e), 
for different values of N and Xu = ln2 (a), lnl.3 (b). Lines 
correspond to analytical predictions; symbols to numerical 
simulations using the Bernoulli (a) or triangular (b) maps. 
Synchronization is transversally stable in the region between 
the couple of curves for each set of values of the parameters. 



that for w = 0.074 yields Xjj ~ lnl.3 (notice that 
Xu = —wlnw — (1 — w)ln[l — w]). Additional tests 
(results not shown here) were performed for other maps 
such as the logistic map f{x) = fix{l — x) with /i = 4 
(hence Xu = In 2) and /x = 3.6533 (hence Xu — In 1.30) 
yielding the same degree of agreement. For these interval 
maps, in principle, one must have < £ < 1 in order to 

(i) 

guarantee that the state variables a;„ will remain inside 
the interval [0, 1]. But, reinjection into the interval can 
be performed through an operation, for instance (mod 
1), such that it does not spoil the Lyapunov exponent 
of the chaotic uncoupled map. If trajectory points were 
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not reinjected, one can still look at our results as valid 
for trajectories or trajectory segments as long as the state 
variables remain confined within the given interval. Any- 
way, for other maps such as f{x) = exp{— [(a; — 0.5)/c7]^} 
[15], one may have any coupling e > since the map is 
naturally defined in the full real axis. Tests performed 
with this Gaussian map (results not shown in this letter) 
are also in good accord with theoretical predictions. 

In general terms, wc observe that for weak coupling, 
the maps do not synchronize. As the coupling strength 
increases, synchronization can occur depending on the 
lattice parameters (a, N, Xjj). However, a too high cou- 
pling intensity e > e'^ has a destabilizing influence on the 
CS state and the lattice no longer synchronizes. 

Concerning lattice size, already Fig. 1 exhibits the in- 
tuitive fact that it is more difficult to synchronize a larger 
lattice than a shorter one, all other parameters being kept 
fixed. In the limit N oo we obtain 



Sc{a, 00, Xu) 



1 



1 - C{a) ' 



(8) 



where C{a) = limjv^oo b'^^^ /rj. This limit is equal to one 
for a > 1, so that Eq. (8) yields a divergent result. For 
a outside the domain of convergence of the series, i.e. 
a <1, one gets 



C{a) 



1 



a 



^ cos(a;) 



dx. 



In that same range of a one has 

e'^{a,oo,Xu) = l-|-e" 



(9) 



(10) 



which is independent on a, thus it yields a straight line in 
the plots of Fig. 1. From the intersection of eda, oo, Xu) 
with £^(a, oo,A!7) it results a critical value of the inter- 
action range ac, such that, for a < ac < 1 (ac < d in 
the d-dimensional case [9]), synchronization is possible 
even in the thermodynamic limit N —> oo for an appro- 
priate window of s. Observe the corresponding domains 
in Fig. 1. 

As N diminishes, the upper curve in Fig. 1, which is 
a straight line for infinite A'^, gains a negative inclination 
and extends for large a to values of e less than unity, indi- 
cating that the destabilizing effect of very strong coupling 
is more easily attained. The lower curve segment has a 
positive inclination, connecting to the upper curve at a 
point that forms a cusp for small N. For still smaller A^ 
(e.g., A^ < 5 for Ai/ = ln2 and TV < 8 for Xu = In 1.3) 
the two critical curves do not intersect each other even 
in the limiting case of first neighbors (a — *■ oo). 

The effect of lattice size can also be observed in Fig. 2 
that exhibits the synchronization domains in the plane 
{N, e) for different values of a. Similar plots have been 
observed in scale- free networks [16], as if there were an 
average or effective a in such cases. While for a > ac 
there is un upper bound Nh{a, Xu) of the number of maps 
for which synchronization occurs; for a < Uc any number 



of maps synchronize (because the critical curves in Fig. 2 
do not intersect). Gencrically it is easier to synchronize 
a small number of maps. Consistently with this obser- 
vation, lattices of small size (e.g., A^ < 5 for A[/ = ln2 ) 
can synchronize for any a, for a certain window of e that 
narrows with increasing a. For A^ < 3 there is, naturally, 
no dependence on a and the lattice syncronizes for any 
Xu > 0. 
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FIG. 2. Synchronization critical lines in the plane (A'^, e) 
for different values of a and Xu ~ In 2. In this case Uc « 0.77. 
Symbols correspond to theoretical calculations. Dotted lines 
are guides to the eyes. Synchronization is transversally stable 
in the region between the two curves for each set of values of 
the parameters. 



Although Figs. 1(a) and 1(b) yield qualitative similar 
results, their comparison makes clear that, as expected, 
the more chaotic the uncoupled maps are, the more diffi- 
cult becomes to obtain their synchronization. The influ- 
ence of the Lyapunov exponent Xu on the synchroniza- 
tion domains in the parameter space {a,e) is displayed 
in Fig. 3. The exhibited numerical results were acquired 
for A^ = 21 either Bernoulli or triangular maps. 

If the positive Lyapunov exponent of the uncoupled 
map increases, the synchronization domain shrinks, col- 
lapsing in the limit Xu — > oo. In the opposite limit of 
Xu 0+ one gets 





2(1 



V 



and 

-1 



(11) 

(12) 



This limit value of e'^ depends on a and A^. If A^ ^ oo and 
a 0(oo), £^ goes to 2.0 (1.0) in the limit of vanishing 
chaos. (Although N = 21, all these extreme behaviors 
are already insinuated in Fig. 3.) As a consequence, the 
critical value ac < 1 below which the lattice synchro- 
nizes in the thermodynamic limit depends on the degree 
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of chaoticity of the uncoupled maps. This dependence 
can be explicitly obtained by inversion of 



\u = In 



C{ac) 



(13) 



extracted from Eqs. (8) and (10). The critical value ac 
as a function of Xu is displayed in Fig. 4. In the limit of 
vanishing (infinite) Xu, olc goes 1.0 (0.0). 




FIG. 3. Synchronization domains in parameter plane (a, e) 
for 21 coupled chaotic maps and various values of Xu- Nu- 
merical simulations were performed for Bernoulli or triangular 
maps. Critical curves were obtained analytically from Eqs. (5) 
and (6). 
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FIG. 4. Critical value Uc, below which synchronization is 

stable even in the thermodynamic limit, as a function of \jj 
(symbols), determined from Eq. (13). The dotted line is a 
guide to the eyes. 



When a > ac, it must be < Ni,{a,Xu) for the lat- 
tice to synchronize, where Nt decreases with increasing 
a — ac (as shown in Fig 2). In the limit a — > oo, it is 
easy to obtain, from the condition e'^ > Sc, an approxi- 
mate expression for the maximal size, valid when one has 
sufficiently small Xu and large Ni,: 



Nmax = Nb{(X), Xu) < IT 



(14) 



In Fig. 5, we exhibit the maximal size N^ax for which 
synchronization can be achieved in the limit of nearest- 
neighbor couplings (hence for any a) as a function of Xu, 
together with the approximation given by Eq. (14). 
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FIG. 5. Maximal lattice size N^ax = Nb{a,Xu), for which 
synchronization can be achieved for any a, as a function of the 
chaoticity indicator \u (circles), determined from the condi- 
tion e'c > Ec for a = oo. The solid line corresponds to the 
approximation given by Eq. (14). 



Summarizing, we have presented numerical and ana- 
lytical results for the CS states in lattices of coupled 

identical chaotic maps with interactions that decay with 
distance as a power law. We have scrutinized the role 
of the system parameters in the ability of the lattice 
to attain complete synchronization using various chaotic 
maps. We observed, in the coupling parameters plane, 
an overall decrease of the area of the synchronization re- 
gions, as the number of coupled maps is increased. More- 
over, the shape of these regions is bounded by critical 
curves which vary with the lattice size in a fashion we 
were able to predict analytically based on the Lyapunov 
spectrum of the synchronized state of the lattice, for vir- 
tually any chaotic map, with excellent agreement with 
numerical results. Similar findings have been given for 
the dependence of the synchronization regions on the de- 
gree of chaoticity of uncoupled maps, as well as for the 
maximal lattice size for which complete synchronization 
of chaos can be achieved. 
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